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Abstract 


We study the possible existence of nonstrange dibaryons and tribaryons by 
solving the bound-state problem of the two- and three-body systems composed 
of nucleons and deltas. The two-body systems are NN, NA, and A A, while 
the three-body systems are NNN, NNA, AAA, and AAA. We use as 
input the nonlocal NN, NA, and AA potentials derived from the chiral 
quark cluster model by means of the resonating group method. We compare 
with previous results obtained from the local version based on the Born- 
Oppenheimer approximation. 
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I. INTRODUCTION 


Systems without strangeness are those which involve only nucleons and nonstrange 
mesons like the pion or the eta. In a series of recent papers the suggestion has been made 
that it may be possible to observe unstable nonstrange two- and three-baryon states cor¬ 
responding to the bound-state solutions of the various systems composed of nucleons and 
deltas 0-§]. These are the systems NN, NA, AA, NNN, NNA, NAA, and AAA. The 
bound states involving one or more unstable particles will show up in nature as dibaryon or 
tribaryon resonances. In the case of two-body systems (dibaryons) they will decay mainly 
into two nucleons and either one or two pions, while for the three-body case (tribaryons) 
they will decay mainly into three nucleons and either one, two, or three pions. 

In the previous calculations of our group [Q-^, the Born-Oppenheimer approximation was 
used in order to obtain a local potential for the baryon-baryon interactions {NN, NA, AA). 
In this work, we will overcome the Born-Oppenheimer approximation by working directly 
with a nonlocal potential derived within the resonating group method (ROM) formalism. 
This method allows, once the Hilbert space for the six-body problem has been hxed, to treat 
the inter-cluster dynamics in an exact way. 

In order to perform the NNN, NNA, AAA, and AAA calculations we will take ad¬ 
vantage of the experience gained in the three-nucleon bound-state problem In that 

case one knows that the dominant conhguration of the system is that in which all particles 
are in S-wave states. However, in order to get reasonable results for the binding energy, the 
S-wave two-body amplitudes used as input in the Faddeev equations must contain already 
the effect of the tensor force. Thus, for example, in the case of the Reid soft-core potential 
if one considers only the S-wave conhgurations but neglects the tensor force in the two-body 
subsystems the triton is unbound. However, if one includes the effect of the tensor force 
in the nucleon-nucleon channel, but uses only the and ^Si components of the 

two-body amplitudes in the three-body equations (2-channel calculation), one gets a triton 
binding energy of 6.58 MeV. Notice that including the remaining conhgurations (34-channel 
calculation), leads to a triton binding energy of 7.35 MeV [^. This means that the S-wave 
truncated T-matrix approximation leads to a binding energy which differs from the exact 
value by less than 1 MeV. Therefore, by means of our approach we will not study exact 
binding energies but which are the best candidates for bound states and the ordering of the 
different NNN, NNA, NAA, and AAA states. 

NN, NA, and AA interactions have been derived in the past in the framework of meson- 
exchange models or phenomenological potentials These models have been used over the 
years to ht the NN data very accurately. However, in the NA and A A sectors experimental 
data are so scarce that it is not possible to obtain reliable values of the parameters involved 
in the interaction. The situation is different in the case of quark cluster models |T0|-[I^. In 
these models the basic interaction is at the level of quarks involving only a quark-quark-held 
(pion or gluon) vertex. Therefore its parameters (coupling constants, cutoff masses, etc.) 
are independent of the baryon to which the quarks are coupled, the difference among them 
being generated by SU(2) scaling, as explained in Ref. [|^. Moreover, quark models provide 
a dehnite framework to treat the short-range part of the interaction. The Pauli principle 
between quarks determines the short-range behavior of the different channels without addi¬ 
tional phenomenological assumptions. In this way, even in the absence of experimental data. 
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one has a complete scheme which starting from the NN sector allows us to make predictions 
in the NA and A A sectors. This fact is even more important if one takes into account that 
the short-range dynamics of the NA and A A systems is to a large extent driven by quark 
Pauli blocking effects, that do not appear in the NN sector. Pauli blocking acts in a selec¬ 
tive way in those channels where the spin-isospin-color degrees of freedom are not enough 
to accommodate all the quarks of the system Therefore, meson-exchange models 

cannot fully include the effect of quark Pauli blocking through its purely phenomenological 
short-range channel-independent part. 

The lifetime of the bound states involving one or more deltas should be similar to that 
of the A in the case of very weakly bound systems and larger if the system is very strongly 
bound. Therefore, these dibaryon and tribaryon resonances will have widths similar or 
smaller than the width of the A so that, in principle, they are experimentally observable. 
Also, we want to emphasize that the possible detection of dibaryon and tribaryon resonances 
does not constitute an exotic subject since, in principle, any nucleus with at least three 
nucleons can serve as test system that may be excited by forming a tribaryon 0. 

The paper is organized as follows. In section II we present the basic quark-quark inter¬ 
action and we describe the method to obtain the resonating group method baryon-baryon 
potentials. Section III is dedicated to discuss the formalism to solve the bound state problem 
for the cases of systems of identical and non-identical particles, respectively. In section IV 
we give our results and we present the conclusions in section V. 


II. THE TWO-BODY INTERACTIONS 


The basic two-body interactions, Vab^ab-, between baryons A and B that are going to be 
used in this work are the nucleon-nucleon interaction V^n^nn, the nucleon-delta interaction 
Vna^na, and the delta-delta interaction Vaa^aa- These baryon-baryon interactions were 
obtained from the chiral quark cluster model developed elsewhere |^. In this model baryons 
are described as clusters of three interacting massive (constituent) quarks, the mass coming 
from the breaking of chiral symmetry. The ingredients of the quark-quark interaction are 
conhnement, one-gluon (OGE), one-pion (OPE) and one-sigma (OSE) exchange terms, and 
whose parameters are hxed from the NN data. Explicitly, the quark-quark [qq) interaction 
is, 

~ Vcon{rij) + VoGE{rij) + VbpE(f^p') + PosE(f^p'), (1) 

where fij is the ij interquark distance and 
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Although taken to be linear for consistency with the baryon and meson spectra, the detailed 
radial structure and strength of the conhning potential is meaningless for the two-baryon 
Oc is the conhnement strength, the A’s are the SU(3) color matrices, the 
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a’s (r’s) are the spin (isospin) Pauli matrices, Sij is the usual tensor operator, (tUtt, ma) 
is the quark (pion, sigma) mass, is the gg-gluon coupling constant, ach is the gg-meson 
coupling constant and A is a cut-off parameter. 

For the present study we make use of the nonlocal potentials derived through a 
Lippmann-Schwinger formulation of the RGM equations in momentum space 0- The 
formulation of the RGM for a system of two baryons, Bi and B 2 , needs the wave function 
of the two-baryon system constructed from the one-baryon wave functions. The two-baryon 
wave function can be written as: 


= A[x{P)4>BAnBMB2{'PiB^)xmB2ic[^\ (V 

where A is the antisymmetrizer of the six-quark system, x(P) is the relative wave-function of 
the two clusters, (pBiip^gJ is the internal spatial wave function of the baryon Bi, are the 
internal coordinates of the three quarks of baryon Bi. X^B 2 denotes the spin-isospin wave 
function of the two-baryon system coupled to total spin (S') and isospin (T), and hnally, 
^c[2^] is the product of two color singlets. 

The dynamics of the system is governed by the Schrodinger equation 


{H - Pt)|^ >= 0 (5T|(7f - Et)\^ >= 0, 


( 8 ) 


where 

N W 2 

(9) 

i=l 

being Tc,m. the center of mass kinetic energy, Vij the quark-quark interaction described above, 
and mq the constituent quark mass. 

Assuming the functional form 

//, 2 \ 3/4 

j , (10) 

where b is related to the size of the nucleon quark core, Eq. (j^) can be written in the 
following way, 

X{P') + J K{P', PMP^)dP^ = 0 ( 11 ) 
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where 


K{P\P^) 


RGM 


Vo(P',Pi) + ^'^VEx(P',Pi) 


( 12 ) 


contains the direct and exchange RGM potentials, the later one coming from quark anti¬ 
symmetry. K{P',Pi) is the nonlocal potential. /,From Eq. o a set of coupled Lippmann- 
Schwinger equations can be obtained and solved using standard techniques. The parameters 
of the model are summarized in Table I. They have been hxed in order to obtain the best 
£t of the two-nucleon sector (deuteron binding energy and S-wave NN scattering phase 
shifts) and the A — N mass difference. In particular, the mass of the quark {mg) is taken 
to be 1/3 of the nucleon mass. The pion mass (m^) is its experimental value. The chiral 
coupling constant (cxc/i) has been determined to reproduce the long-range OPE interaction 

and is given by ach = , where the ttNN coupling constant is taken to be 


2 

9tvNN 


=13.87. The sigma mass is hxed by the chiral symmetry relation ~ -|- (2mg)^. 

The parameter b, which determines the size of the nucleon quark content, was determined 
by comparing the adiabatic NN potential calculated from the wave function solution of the 
bound state problem for the potential given by Eq. (1) to the NN potential calculated with 
a single gaussian of parameter b. A, which controls the pion-gluon proportion in the model 
and, as a consequence, the strength of the tensor force, has been taken to reproduce the 
deuteron binding energy in the presence of AA channels. As the OPE provides part of the 
A — N mass difference, the value of the strong coupling constant (a^) is determined to ob¬ 
tain the remaining A — N mass difference. Finally, the value of Oc, quoted for completeness 
because its contribution to the baryon-baryon potential is negligible 
the stability condition for the nucleon = 0 |TT[ . 
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III. INTEGRAL EQUATIONS 

We will describe in this section the formalism required in the cases of the two-body 
systems NN, NA, and AA and the three-body systems NNN, NNA, NAA and AAA. 


A. The two-body systems 


If we consider two baryons A and R in a relative S-state interacting through a potential 
that contains a tensor force, then there is a coupling to the AB D-wave so that the Lippmann- 
Schwinger equation of the system is of the form 


‘ {PuPi] E) = ' (PuPi) + Pi dp'i VPP; ^{Pi.p'i) 

1 


rhsil'i s'f 




X- 




E - pf/2r]i + ie 


(13) 


where ji and A are the angular momentum and isospin of the system, while ksi, l[s[, and I'-s" 
are the initial, intermediate, and hnal orbital angular momentum and spin of the system, 
respectively. Pi and pi are, respectively, the relative momentum and reduced mass of the 
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two-body system. We give in Tables II, III and IV the corresponding NN, NA and AA 
two-body channels in a relative S-wave that are coupled together for the two possible values 
of j and i (since the NN state is the one with the lowest mass, in the case of this system we 
have considered also the possibility of transitions to higher mass states like NA and AA). 
In the cases of the NN and AA systems the Pauli principle requires that (—= —i. 

As mentioned before, for the solution of the three-body system we will use only the 
component of the T-matrix obtained from the solution of Eq. (^) with li = /" = 0, so that 
for that purpose we dehne the S-wave truncated amplitude 




( 14 ) 


B. The three-body systems 


If we restrict ourselves to the conhgurations where all three particles are in S-wave states, 
the Faddeev equations for the bound-state problem in the case of three particles with total 
spin S and total isospin / are 


Sjij 


q^jdqj / dcoseti-siiiiPhPi] E - q}/2ui) 


1-1 
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X- 
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( 16 ) 


where pi and g* are the usual Jacobi coordinates and pi and z/j the corresponding reduced 
masses 
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with ijk an even permutation of 123. The momenta p[ and pj in Eq. (P!3D are given by 
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are the spin-isospin coefficients 
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where W is the Racah coefficient and cxj, Sj, and S {ji, ii, and I) are the spin (isospin) of 
particle i, of the pair jk, and of the three-body system, respectively. 
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Since the variable pi, in Eqs. (|^ and (^), runs from 0 to cx), it is convenient to make 
the transformation 


Xi = 


Pi - d 
Pi + d' 


( 21 ) 


where the new variable Xi runs from —1 to 1, and d is a scale parameter. With this trans¬ 
formation Eq. (p!5D takes the form 
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Since in the amplitude ti-siuixi^x'p^e) the variables Xi and x'i run from —1 to 1, one can 
expand this amplitude in terms of Legendre polynomials as 




(23) 


where the expansion coefficients are given by 




2n -|- 1 2 m -|- 1 
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Applying expansion (^) in Eq. (^) one gets 
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where satishes the one-dimensional integral equation 
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The three amplitudes T[^gj{qi), T^sf^{q 2 )i and T^^g^^^iqs) in Eq. (^ are coupled to¬ 
gether. The number of coupled equations can be reduced, however, since some of the particles 
are identical. In the case of three identical particles [NNN and AAA systems) we have 
that all three amplitudes are equal and therefore Eq. (p6D becomes in this case. 


TSr(qi) = 2E f <'9j2i”3r’‘'(<6,®;-B)r™’‘’(®). 


(28) 
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We give in Table V the three NNN states characterized by total spin and isospin (S', I) that 
are possible as well as the two-body NN channels that contribute to each state. In Table VI 
we give the 25 AAA states characterized by total spin and isospin (S', I) that are possible 
as well as the two-body AA channels that contribute to each state. 

In the case where two particles are identical and one different {NNA and NAA systems) 
two of the amplitudes are equal. The reduction procedure for the case where one has two 
identical fermions has been described before and will not be repeated here. With the 

assumption that particle 1 is the different one and particles 2 and 3 are the two identical, 
only the amplitudes and T^|*^(g 2 ) are independent from each other and they 

satisfy the coupled integral equations 

r;s;‘(9i) = 2E /“rf93A5sir“(9i.93;B)rs»(«s), ( 29 ) 

713212 ^ 
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with the identical-particles phase 

Iden = I + ai + a'i - S 2 + Ti + T'i - i 2 . (31) 

Substitution of Eq. (^) into Eq. (^) yields an equation with only the amplitude T 2 

_ POO 
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We give in Table VII the 9 NNA states characterized by total spin and isospin (S', I) 
that are possible as well as the two-body NA and NN channels that contribute to each 
state. In Table VIII we give the 16 AAA states characterized by total spin and isospin 
(S', I) that are possible as well as the two-body NA and AA channels that contribute to 
each state. 


C. Numerical solutions 


In order to hnd the bound-state solutions of Eqs. (0. @ and (|3^ we drop the 
inhomogeneous term in Eq. ( 0 ) (of course, in the solution of the three-body problem we 
use as input the solutions of the inhomogeneous Eq. (pA])) and replace the integral by a 
sum applying a numerical integration quadrature . In this way, Eqs. ( 0.0 and (|32[) 
become a set of homogeneous linear equations. This set of linear equations has solutions 



only if the determinant of the matrix of the coefficients (the Fredholm determinant) vanishes 
for certain energies. Thus, the procedure to hnd the bound states of the system consists 
simply in searching for the zeroes of the Fredholm determinant as a function of energy. We 
checked our program by comparing with known results for the three-nucleon bound-state 
problem with the Reid soft-core potential 0. We found very stable results taking for the 
scale parameter d = 3 fm“^, a number of Legendre polynomials L = 10, and a number of 
Gauss-Legendre points N = 12. 


IV. RESULTS 


We will now present the results of our nonlocal calculations for the 7 systems correspond¬ 
ing to the two- and three-body bound-state problem of nucleons and deltas, and compare 
them to previous calculations which have been done by our group based on the local poten¬ 
tials obtained from the Born-Oppenheimer approximation. 

The two body interaction in the NA states {j,i) = (1,1) and (2,2), and those of the 
A A states (j, i) = (2,3) and (3,2) present quark Pauli blocking. As a consequence, a strong 
repulsive core appears in the baryon-baryon potential. The reason for that is based on the 
fast decrease of the norm of the six-quark wave function when R —>■ 0 . A similar analysis 

performed in terms of the SU(4) symmetry shows the presence of a forbidden state. From 
the physical point of view, it is connected with the lack of enough degrees of freedom to 
accommodate all the quarks. It is important to note that the origin of this repulsion is not 
the same as in the NN channels, because they do not show a forbidden state but a mixing 
of [6] with the [4, 2] six-quark orbital symmetry. Technically, the reason for such a strong 
repulsive core is the presence of nodes in the inner region of the relative wave function of Eq. 
(^). This behavior originates essentially from the condition that the relative wave function 
should be orthogonal to the forbidden state due to the Pauli principle |^. The forbidden 
state should then be eliminated from the relative wave function for each partial wave. This 
procedure is tedious both from the conceptual and numerical point of view |^J^ . It has 
been demonstrated that for the Pauli blocked channels the local NA and AA potentials 


reproduce the qualitative behavior of the RGM kernels after the subtraction of the forbidden 
states. This is why we used in our calculations the local version of the quark Pauli blocked 
channels mentioned above. 

In the case of the three-body systems we calculated the binding-energy spectrum (that 
is, the energy of the states measured with respect to the three-body threshold) as well as 
the separation-energy spectrum (that is, the energy of the states measured with respect to 
the threshold of one free particle and a bound state of the other two). The deepest bound 
three-body state is not the one with the largest binding energy but the one with the largest 
separation energy, since that state is the one that requires more energy in order to become 
unbound (that is, to move it from the bound state to the nearest threshold). 


A. The NN system 

We found that of the two states of Table II only the one with (j, i) = (1, 0), that is the 
deuteron, is bound. The nonlocal model gives a deuteron binding energy of 2.14 MeV, while 
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the local version gave an energy of 3.13 MeV. These results are shown in Table IX. 

The exact chiral quark cluster model NN potential |]^] gives a deuteron binding energy of 
2.225 MeV. This value was obtained by taking into account the AA partial wave {Iaa, saa ) = 
(4, 3) coupled together in addition to those given in Table II. Since in our calculation we 
consider only S- and D-waves, we omit the AA {Iaa, saa) = (4, 3) partial-wave contribution, 
and we obtain instead a deuteron binding energy of 2.14 MeV, which differs less than 0.1 
MeV from the exact calculation. 


B. The NA system 


We give in Table X the results for the binding energies of the NA system. Out of the 
four possible NA states of Table III only one, the {j,i) = (2,1), has a bound state which 
lies exactly at the NA threshold for the local model. However, if we use the nonlocal 
model we hnd instead a bound state of 0.141 MeV. The states {j,i) = (1,1) and (2,2) are 
unbound because they present quark Pauli blocking [|T^ and therefore they have a strong 
repulsive barrier at short distances in the S-wave central interaction. These two states play 
an important role in the three-body spectrum. The state {j,i) = (2,1) can also exist in 
the NN system and there it corresponds to the ^D 2 partial wave which has a resonance at 
an invariant mass of 2.17 GeV [p^ -p6|. This means that the NA bound state may decay 
into two nucleons and appear in the NN system as a resonance. The NA bound state has 
for both local and nonlocal models energies very close to the NA threshold, so that the 
invariant mass of the system is also very close to 2.17 GeV. Thus, one or another of our 
models predict the NN ^D 2 resonance as being a NA bound state. 


C. The A A system 

We give in Table XI our results for the AA system. Out of the eight possible AA states 
given in Table IV with nonlocal interactions hve have a bound state, whereas the local 
interactions bind six of them (in both local and nonlocal models there are no excited states in 
any of the channels). It is interesting to note that the predicted bound states: (j, i) = (1, 0), 
(0,1), (2,1) and (3,0), also appear in the case of the NN system. In the nonlocal model, 
we hnd that the deepest bound state is (j,z) = (1,0), the second {j,i) = (0,1), the third 
(jG) = (3,0) and the fourth {j,i) = (2,1). This clearly shows that there is a qualitative 
similarity between the AA and NN systems (both are systems of identical particles). Three 
of these states appear also in the case of the NN system. The (j, z) = (1,0) state is of 
course the deuteron, the (j, z) = (0,1) is the virtual state and the (j, z) = (2,1) state 
is the ^D 2 resonance that lies at ~ 2.17 GeV ||^ (note that the NN resonance has 
no counterpart in Table XI because we calculated only even-parity states and has odd 
parity). Thus, the (j, z) = (3,0) state, which is also allowed in the case of the NN system, 
would correspond to a new nucleon-nucleon resonance that is predicted by our model. The 
(i, G = (3,0) channel corresponds in the case of the NN system to the partial wave. 
Some indication of the (3, 0) resonance can already be seen in the most recent analysis of 
the NN data by Arndt et al. PB| . 
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The channels (j, i) = (2, 3) and (3, 2) are nnbonnd becanse they have a strong repnlsive 
barrier at short distances in the S-wave central interaction. This strong repnlsion originates 
from the qnark Panli blocking prodnced by the satnration of states that occnrs when the 
total spin and isospin are near their maximnm valnes |^. As we will see later in the 
discnssion of the AAA resnlts, these repnlsive cores in the (3, 2) and (2, 3) channels largely 
determine the three-body spectrnm. 

^From Table XI we note that the two-body AA bonnd states which have low qnantnm 
nnmbers are deeper for the nonlocal model than with the local one. This pecnliar featnre 
resnlts to be conversely for the case of high qnantnm nnmbers. 


D. The NNN system 

As another test of the reliability of onr model in the case of the three-baryon system 
we solved the NNN bonnd-state problem. We fonnd that of the states of Table V only 
the state with {S,I) = (|,|), that is the triton, has a bonnd state. By nsing the local 
potentials we obtain a binding energy of 5.76 MeV for the triton. On other hand, if we nse 
the nonlocal potentials as inpnt we find a triton binding energy of 6.52 MeV. For comparison, 
we notice that the triton binding energy for the Reid-soft-core potential in the truncated 
T-matrix approximation is 6.58 MeV. Since the experimental value is Bexp = 8.49 MeV 
the difference with our theoretical result, of about 3 MeV, is a measure of the uncertainty 
of our calculation in the case of the three-baryon system. We show in Table XII the results 
of our calculations for the NNN system. There, is the binding energy of the system 
and i ?3 — i ?2 is the separation energy, being B 2 the binding energy of the deepest bound 
two-body channel that contributes to the three-body state (see Table IX). 


E. The NNA system 


We show in Table XIII the results of our calculations for the NNA system. 

One may have hoped to hnd several bound states in this system, due to the fact that 
the NA two-body subsystem has a bound state in the channel (j, i) = (2,1) and the NN 
two-body subsystem has a bound state in the channel (j, i) = (1,0) and an almost-bound 
state in the channel {j,i) = (0, !)• This is not the case however, and as a matter of fact, 
with the nonlocal potentials as input only two of the 9 possible three-body states given 
in Table VII are bound. Because of the attractive contribution of the NA {j,i) = (2,1) 
bound state with the nonlocal model, the three-body state (|, |) results to be very weakly 
bound, at an energy of 0.143 MeV, and a separation energy scarcely different from zero. 
That means that the (5,/) = (|, |) state is very near the NNA threshold and therefore 
it represents the tribaryon resonance with the lowest possible mass since it can decay into 
three nucleons and one pion. Also, for this case the three-body state (|, |) is bound. As it 
can be seen from Table VII, this state has the contribution of all the two-body NA and NN 
channels. In spite of the fact that the NA two-body channels (j, i) = (1,1) and (2, 2) present 
Pauli blocking |14], and therefore they have a strong repulsive barrier at short distances in 


the S-wave central interaction, the attractive contribution of the NA {j,i) = (2,1) and NN 
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(j) i) = (1) 0) channels results to be enough to weakly bound this state. We note that neither 
one of the three-body states (S', I) = (|, |) and (|, |) is bound with local interactions. 

F. The A^AA system 

The results for the NAA system are shown in Table XIV. Similarly to the case just 
discussed, in our calculations with nonlocal interactions we found that three of the 16 possible 
7VAA states given in Table VIII are bound. They are the (S',/) = (|, |),(|, |), and (|, |) 
states and their corresponding bound state energies are 0.630 MeV, 8.158 MeV, and 0.181 
MeV, respectively. In the case of the states (S',/) = (i, |) and (S',/) = (|, \) the repulsive 
barrier due the quark Pauli blocking in the NA states (j, i) = (1,1) and (2, 2) is less strong 
than the attraction due to the state (j, i) = (2,1), so that they result to be bound states in 
the nonlocal model. The state (S', /) = (|, |) is the weakest bound state of this system, since 
in addition to the contribution of the NA quark Pauli blocking channels, there exists that of 
the A A quark Pauli blocking channels (j, i) = (2, 3) and (3, 2). This conhrms what we have 
mentioned before that it is the structure of the interaction of the two-body system the one 
which largely determines the three-body spectrum. Thus, the nonlocal interactions predict 
the bound states (S',/) = (|, |), (|, |) and (|, |), which in principle may be observable as 
tribaryon resonances which decay into three nucleons and two pions with masses close to 
the NAA threshold. 


G. The AAA system 

We show in Table XV the results for the AAA system. The system has 4 bound states 
while by using the local interactions the system had 7 bound states. /^From Table XV we 
observe that the three states which are missing in the nonlocal version are barely bound in 
the local version, i.e., they have very small separation energies. Since the nonlocal interaction 
tends to lower the attraction in all the AAA channels it is not surprising that those which 
were barely bound have now disappeared. The more strongly bound three-body state (that 
is, the one with the largest separation energy) is the (S', /) = (|, |) state which has precisely 
the quantum numbers of the triton. This shows again, like in the AA and NN systems, the 
similarity between the AAA and NNN systems. 

The reason why the (S',/) = (|, |) state is the more strongly bound is very simple. As 
shown in Table VI, this is the only state where none of the two-body channels with a strong 
repulsive core (j, i) = (2, 3) or (3, 2) contribute. In all the other three-body states the strong 
repulsion of the (j, i) = (2, 3) and (3, 2) channels either completely destroys the bound state 
or allows just a barely bound one. The state (S',/) = (|, |) comes next with respect to 
separation energy. This state (S',/) = (|, |) has a somewhat anomalous behavior since it 
has a relatively large separation energy. This behavior is sort of accidental and it can be 
understood as follows. As seen in Table VI, there are 4 two-body channels contributing to 
the (S', /) = (|, |) state, the two attractive ones (j, i) = (2,1) and (3, 0) and the two repulsive 
ones (j, i) = (2,3) and (3,2). However, as one can see in Table XI the attractive channels 
(2,1) and (3, 0) have bound states dX E = —7.4 MeV and E = —7.8 MeV, respectively, 
for the nonlocal version, and E = —30.5 MeV and E = —29.9 MeV, respectively, for the 
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local version, so that the poles in the scattering amplitudes of these two channels are very 
close together and therefore there is a reinforcement between them, which gives rise to the 
anomalously large separation energy in both versions. 


V. CONCLUSIONS 


By using both the local and nonlocal models we have studied the bound-state solutions 
of the two- and three-body systems composed of nucleons and deltas. First of all we would 
like to emphasize the goodness of the Born-Oppenheimer approximation, producing results 
very similar to the usually more involved RGM results. We conclude that the more realistic 
nonlocal interactions produce in the two-body systems NN, NA, and AA one, one, and hve 
bound states respectively. The bound states of the unstable systems NA and A A correspond 
to dibaryon resonances that decay mainly into two nucleons and one pion and two nucleons 
and two pions, respectively. The NA bound state with {j,i) = (2,1) and M Ri 2.17 GeV is 
the dibaryon resonance with the lowest possible mass and the one which seems to be well 
conhrmed by experiment. The hve A A bound states of the nonlocal potentials correspond 
to dibaryon resonances with masses between 2.4 and 2.5 GeV. The (j, i) = (3,0) AA state 
would correspond to a new nucleon-nucleon resonance predicted by our model. A possible 
signal of this resonance appears in a recent analysis of NN data up to 3 GeV by Arndt et 
al. 1^. With respect to the three-body systems we found that the NNN has one bound 


state, the AAA has four bound states, the NN A has two bound states, and the N AA 
has three bound states. The predicted NN A states with (S', I) = (|, and (S, I) = (|, |) 
which correspond to M k, 3.4 GeV are the tribaryon resonances with the lowest mass and 
therefore the ones that would be more easy to detect experimentally. 
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TABLES 


TABLE 1. Quark model parameters. 


mq(MeV) 

313 

6(fm) 

0.518 


CXg 

0.498 

(MeV-fm-1) 

67.0 

^ch 

0.027 

m(j(fm“^) 

3.513 

m^(fm“^) 

0.70 

A(fm“^) 

4.3 


TABLE II. NN channels {Inn, snn), A'A channels {Ina, sna), and AA channels {Iaa,s/\a) 
that are coupled together in the and ^S'o NN states. 

NN state _j_ i _ {Inn, snn) _ {Ina, sna) _ {Iaa, saa) 

1 0 (0,1),(2,1) - (0,1),(2,1),(2,3) 

0 1 ( 0 , 0 ) ( 2 , 2 ) 


TABLE III. Coupled channels {l,s) that contribute to a given NA state with total angular 


momentum j and isospin i. 


j 

i 

{l,s) 

1 

1 

(0,1),(2,1),(2,2) 

1 

2 

(0,1),(2,1),(2,2) 

2 

1 

(0,2),(2,1),(2,2) 

2 

2 

(0,2),(2,1),(2,2) 


TABLE IV. Coupled channels {l,s) that contribute to a given AA state with total angular 
momentum j and isospin i. 

j 

i 

{l,s) 

0 

I 

(0,0),(2,2) 

0 

3 

(0,0),(2,2) 

1 

0 

(0,1),(2,1),(2,3) 

1 

2 

(0,1),(2,1),(2,3) 

2 

I 

(0,2),(2,0),(2,2) 

2 

3 

(0,2),(2,0),(2,2) 

3 

0 

(0,3),(2,1),(2,3) 

3 

2 

(0,3),(2,1),(2,3) 
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TABLE V. Two-body NN channels (j, i) that contribute to a given NNN state with total 


spin S and isospin I. 


5 

I 

0/ *) 

1/2 

1/2 

(1,0),(0,1) 

1/2 

3/2 

(0,1) 

3/2 

1/2 

(1,0) 

TABLE VI. Two-body AA channels (j, i) i 
spin S and isospin I. 

that contribute to a given AAA state with total 

S 

/ 

(j, 0 

1/2 

1/2 

(1,2),(2,1) 

1/2 

3/2 

(1,0),(1,2),(2,1),(2,3) 

1/2 

5/2 

(1,2),(2,1),(2,3) 

1/2 

7/2 

(1,2),(2,3) 

1/2 

9/2 

(2,3) 

3/2 

1/2 

(0,1),(1,2),(2,1),(3,2) 

3/2 

3/2 

(0,1),(0,3),(1,0),(1,2), 

(2,1),(2,3),(3,0),(3,2) 

3/2 

5/2 

(0,1),(0,3),(1,2),(2,1), 

(2,3),(3,2) 

312 

7/2 

(0,3),(1,2),(2,3),(3,2) 

3/2 

9/2 

(0,3),(2,3) 

5/2 

1/2 

(1,2),(2,1),(3,2) 

5/2 

3/2 

(1,0),(1,2),(2,1),(2,3), 

(3,0),(3,2) 

5/2 

5/2 

(1,2),(2,1),(2,3),(3,2) 

5/2 

7/2 

(1,2),(2,3),(3,2) 

5/2 

9/2 

(2,3) 

7/2 

1/2 

(2,1),(3,2) 

7/2 

3/2 

(2,1),(2,3),(3,0),(3,2) 

7/2 

5/2 

(2,1),(2,3),(3,2) 

7/2 

7/2 

(2,3),(3,2) 

7/2 

9/2 

(2,3) 

9/2 

1/2 

(3,2) 

9/2 

3/2 

(3,0),(3,2) 

9/2 

5/2 

(3,2) 

9/2 

7/2 

(3,2) 

9/2 

9/2 
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TABLE VIE Two-body N/S. channels {jNA,iNA) and two-body NN channels {Jnn 
contribute to a given NNA state with total spin S and isospin I. 

, *AfTV) that 

S 

I 

{jNA,iNA) 

ijNN,'iNN) 

1/2 

1/2 

(1,1) 


1/2 

3/2 

(1,1),(1,2) 

(1,0) 

1/2 

5/2 

(1,2) 


3/2 

1/2 

(1,1),(2,1) 

(0,1) 

3/2 

3/2 

(1,1),(1,2),(2,1),(2,2) 

(1,0),(0,1) 

3/2 

5/2 

(1,2),(2,2) 

(0,1) 

5/2 

1/2 

(2,1) 


5/2 

3/2 

(2,1),(2,2) 

(1,0) 

5/2 

5/2 

(2,2) 



TABLE VIII. Two-body NA channels {jNA,iNA) and two-body AA channels ( jaa ,^ aa ) that 
contribute to a given AAA state with total spin S and isospin I. 

S 

I 

{jNA,iNA) 

( j ' aa ,^ aa ) 

1/2 

1/2 

(1,1),(1,2),(2,1),(2,2) 

(1,0),(0,1) 

1/2 

3/2 

(1,1),(1,2),(2,1),(2,2) 

(0,1),(1,2) 

1/2 

5/2 

(1,1),(1,2),(2,1),(2,2) 

(0,3),(1,2) 

1/2 

7/2 

(1,2),(2,2) 

(0,3) 

3/2 

1/2 

(1,1),(1,2),(2,1),(2,2) 

(1,0),(2,1) 

3/2 

2,12 

(1,1),(1,2),(2,1),(2,2) 

(1,2),(2,1) 

3/2 

5/2 

(1,1),(1,2),(2,1),(2,2) 

(1,2),(2,3) 

3/2 

7/2 

(1,2),(2,2) 

(2,3) 

5/2 

1/2 

(1,1),(1,2),(2,1),(2,2) 

(2,1),(3,0) 

5/2 

3/2 

(1,1),(1,2),(2,1),(2,2) 

(2,1),(3,2) 

5/2 

5/2 

(1,1),(1,2),(2,1),(2,2) 

(2,3),(3,2) 

5/2 

7/2 

(1,2),(2,2) 

(2,3) 

7/2 

1/2 

(2,1),(2,2) 

(3,0) 

7/2 

3/2 

(2,1),(2,2) 

(3,2) 

7/2 

5/2 

(2,1),(2,2) 

(3,2) 

7/2 

7/2 

(2,2) 



TABLE IX. Binding energies B 2 of the NN states with total angular momentum j and isospin 
i. B 2 are the results of the local model and B^^ are the results of the nonlocal model. 


j 

i 

Bl'(MeV) 

.B/^^(MeV) 

1 

0 

3.13 

2.14 

0 

1 

unbound 

unbound 


17 












TABLE X. Binding energies B 2 of the NA. states with total angular momentum j and isospin 
i. B 2 are the results of the local model and B^^ are the results of the nonlocal model. 


j 

i 

B^{MeY) 

.B^^(MeV) 

1 

1 

unbound 

unbound 

1 

2 

unbound 

unbound 

2 

1 

0.0 

0.141 

2 

2 

unbound 

unbound 


TABLE XL Binding energies B 2 of the AA states with total angular momentum j and isospin 
i. B 2 are the results of the local model and B^^ are the results of the nonlocal model. 

j 

i 

Bl'(MeV) 

B^^{MeY) 

0 

1 

108.4 

159.5 

0 

3 

0.4 

0.2 

1 

0 

138.5 

190.3 

1 

2 

5.7 

unbound 

2 

1 

30.5 

7.4 

2 

3 

unbound 

unbound 

3 

0 

29.9 

7.8 

3 

2 

unbound 

unbound 


TABLE XII. Binding energies B 3 and separation energies B 3 — B 2 of the NNN states with 
total spin S and isospin I. B 2 and B^ are the results of the local model while B^^ and B^^ are 
the results of the nonlocal model. 


S 

I 

.Bj^(MeV) 

^3 - 

B|'(MeV) 

B|^^(MeV) 

B^^ - B^^iUeV) 

1/2 

1/2 

5.76 


2.63 

6.52 

4.38 

1/2 

3/2 

unbound 


— 

unbound 

— 

3/2 

1/2 

unbound 


- 

unbound 

- 

TABLE XIII. 

Binding energies B^ and separation energies B 3 — 

B 2 of the NNA states with 

total spin S and isospin I. B 2 and B^ are the results of the local model while B^^ and B^^ are 
the results of the nonlocal model. 

S 

I 

B^{MeY) 

1 

cq 

Bl'(MeV) 

B|^^(MeV) 

B^^ - Bf^(MeV) 

3/2 

1/2 

unbound 


- 

0.143 

0.002 

3/2 

3/2 

unbound 


- 

2.280 

0.144 
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TABLE XIV. Binding energies and separation energies — B 2 of the VAA states with 
total spin S and isospin I. and B^ are the results of the local model while B^^ and B^^ are 
the results of the nonlocal model. 


s 

I 

B|'(MeV) 

B^ - Bl'(MeV) 

B|^^(MeV) 

- B^^(MeV) 

1/2 

5/2 

unbound 

- 

0.630 

0.43 

5/2 

1/2 

unbound 

- 

8.158 

0.358 

5/2 

5/2 

unbound 

- 

0.181 

0.04 


TABLE XV. Binding energies B 3 and separation energies B 3 — B 2 of the AAA states with 
total spin S and isospin I. B 2 and B^ are the results of the local model while B^^ and B^^ are 
the results of the nonlocal model. 

S 

/ 

Bj^(MeV) B^ - 

B|'(MeV) 

B^^(MeV) 

B^^ - Bf^(MeV) 

1/2 

1/2 

84.0 

53.5 

16.6 

9.2 

1/2 

3/2 

139.2 

0.7 

unbound 

- 

1/2 

7/2 

6.3 

0.6 

unbound 

- 

3/2 

1/2 

109.5 

1.1 

unbound 

- 

5/2 

1/2 

39.1 

8.6 

9.3 

1.9 

7/2 

1/2 

31.7 

1.2 

7.8 

0.4 

7/2 

3/2 

35.1 

4.6 

9.8 

2.0 
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